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Abstract 



A method of entanglement production is suggested, based on the resonant genera- 

i-rt ! tion of topological modes in systems with Bose-Einstein condensates trapped in optical 

SlJ ' or magnetic lattices. The method makes it possible to regulate the strength of entan- 

^ ■ glement production as well as to govern its time variation. This method can serve as 

a practical tool for quantum information processing and quantum computing. 
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1 Introduction 

Systems with Bose-Einstein condensate (BEC) possess a number of unusual properties (see 
book [1] and review articles [2-7]), which can be employed in a variety of applications. 
One of such important applications is the possibility of creating massive entanglement in 
Bose-condensed systems trapped in lattice potentials [8,9]. The most known examples of 
the lattices are the optical lattices formed by a combination of laser beams [10-12], though 
magnetic fields could also be used for creating magnetic lattices. In experiments, one can 
form optical lattices of varying spacing, depth, and different filling factors, ranging between 
one and 10^ atoms in each lattice site [13,14]. When a lattice is sufficiently deep and the 
filling factor is large, then a lattice site represents a microtrap, in which at low temperature 
BEC can be formed. 

A BEC in a trap possesses a set of discrete or quasidiscrete atomic energy levels. If 
different trapping sites are completely separated from each other, the energy spectrum is 
purely discrete. When atoms can tunnel between the sites, then a line of the discrete 
spectrum widens into a band. In what follows, we consider the situation, when the linewidths 
are not too large, so that the bands are always well separated from each other, that is, when 
the linewidths are much smaller than the spectrum gaps. The latter situation corresponds 
to a quasidiscrete spectrum. Note that we are talking here about the atomic BEC energy 
levels in a trap, which should not be confused with the spectrum of elementary excitations. 

In an equilibrium system, BEC sets in the lowest energy level. But if the system is 
subject to an alternating modulating field, with a frequency in resonance with one of the 
transition frequencies, then a nonground-state BEC can be realized, as was first proposed 
in Ref. [15]. The condensate functions, describing the standard ground-state BEC and the 
nonground-state condensates have different spatial shapes, because of which the condensate 
wave functions, pertaining to different energy levels, can be termed topological modes. The 
properties of these modes have been theoretically investigated in a series of papers [15-35] 
and a dipole topological mode was generated in experiment [36]. A simple example of such 
a mode is a vortex that can be excited in a rotating BEC. 

Dynamics of BEC are usually considered in the frame of the Gross-Pitaevskii equation 
(see Refs. [1,2]), which presupposes the case of zero temperature and very weak atomic 
interactions. The Gross-Pitaevskii equation was also the basis for describing the resonant 
excitation of topological modes [15]. In the present paper, we show that topological modes 
can be created as well at finite temperature, which is due to the resonant mechanism of their 
generation. 

Another our goal is to demonstrate how the regulated generation of topological modes 
in optical lattices can be used for the controlled entanglement production. The possibility of 
effectively varying the state of a complex quantum system and of controlling its entanglement 
are the two key points for realizing quantum information processing and quantum computing 
[37-42]. 



2 Topological modes 

In order to give a general correct definition of topological coherent modes, we need, first, 
to write down the exact equation for the condensate wave function. The latter, for a Bose- 
condensed system, is introduced by means of the Bogolubov shift [43-45] for the field operator 

ij{r,t) -^ ij{r,t)=r]{r,t)+M^,t), (1) 

in which f]{r,t) is the condensate wave function and ■?/'i(r, t) is the field operator of uncon- 
densed particles. The Bogolubov shift (1) explicitly breaks the gauge symmetry of the Bose 
system. It is worth emphasizing that the gauge symmetry breaking is the necessary and suf- 
ficient condition for the occurrence of BEC [46]. After introducing the Bogolubov shift (1), it 
is necessary to resort to a representative statistical ensemble for the Bose system with broken 
gauge symmetry [47,48], defining the appropriate grand Hamiltonian H[iIj] = H[ri,ilJi]. The 
theory of Bose-condensed systems is self-consistent only with the correctly defined grand 
Hamiltonian [49-51]. The general equation for the condensate wave function has the form 

. ^^ = < ^1^ > . (2) 

Assuming the standard energy Hamiltonian with the local interaction potential 

<I>(r) = %b(j) , $0 = 47r ^ , (3) 

vn 

where a^ is the scattering length, we obtain from Eq. (2) the exact equation for the conden- 
sate wave function 

+$0 { [ Po(r, t) + 2pi(r, t) ] ^^(r, t) + ai(r, t)77*(r, t) + ^(r, t) } . (4) 

Here U = f/(r, t) is an external potential and the notation is used for the condensate density 

poir,t) = \r]iv,t)\\ (5) 

the density of uncondensed atoms 

pi(r,t) =<^l(r,t)^i(r,t)>, (6) 

the anomalous average 

cTi(r,t) =<^i(r,t)^i(r,t) >, (7) 

and the triple anomalous average 

e(r,t) =<ijl{r,t)M^,t)M^,t)> . (8) 

Looking for the stationary solutions of Eq. (4) in the common form 

r7„(r,t) = r7„(r)e-'""*. 



we come to the stationary equation for the condensate wave function 
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Vn{r) + 



+$o{ [|77„(r)|2 + 2pi(r)]77„(r) + CTi(r)77„(r)+e(r) } = EnVn{r) 



in which 



Er, 



HO + Un . 



(9) 
(10) 



In equihbrium, BEC corresponds to the lowest energy level 

Eq = min En = Uo , ^o = ^^ u;„ = . 

n n 

But, generally, the nonlinear eigenvalue problem (9) possesses a set of energy levels En and 
of the related eigenf unctions i]n{r), labelled by a multi-index n. The solutions i]n{r) to the 
eigenvalue problem (9) are the topological coherent modes. It is only in the limiting case of 
zero temperature and asymptotically weak interactions, when we can neglect pi(r), (Ti(r), 
and ^(r) in Eq. (9), we come to the stationary Gross-Pitaevskii equation 



— + [/(r) + $ok(r)p 



r]n{r) = En r]n{r) , 
where, actually, solely the ground-state level Eq = /io is to be considered [1]. 



3 Resonant generation 

Equation (9) defines the topological modes as stationary solutions. In order to realize tran- 
sitions between modes, it is necessary to include a time- dependent external potential and to 
consider the temporal equation (4). 

Suppose that at the initial time t = the system is in equilibrium and the condensate 
wave function corresponds to the standard ground-state condensate. 



r]{r,0) = r]o{r) = r]{r) . 
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Let us wish to generate a topological mode labelled by the index n = ni, with the related 
energy Ei = E^- Hence, the transition frequency is given by 



(^10 = El — Eq = En^ — /Uo . 

To generate this mode, it is necessary to apply an alternating field 

V{v, t) = \/i(r) cos(cjt) + 14 (r) sin(cjt) 



(12) 



(13) 



with a frequency u tuned close to the transition frequency (12), so that the resonance 

condition 

A u; 

<1 (A cu = cu - tuio) (14) 



UJ 



be valid. Then the total external potential in Eq. (4) is the sum 

U{r,t) = U{r) + V{r,t) 



(15) 



of a trapping potential and of the alternating field (13). 

The condensate wave function is normalized to the total number of condensed atoms 



No = J\r]{r,t)\^dr = J\r]{r)\^dr. 
Let us introduce a function v'„(r), defined by the equality 



which is normalized to one 



J \^n{r)\^ dr = 1 . 
We shall look for the solution of Eq. (4) in the form 

v{T^,t) =J2 Cn{t) r7„(r) e 



— iuJnt 



(16) 

(17) 
(18) 

(19) 



where u;„ ^ En — Eq is in agreement with Eq. (10). The coefficient function C„(t) is treated 
as a slow function of time, such that 



1 

UJn 



GLLyfi 



dt 



<1. 



(20) 



This allows us to consider C„ as a quasi-integral of motion and to employ the averaging 
method [52] and the scale separation approach [53,54]. For instance, substituting form (19) 
into the normalization condition (16), and averaging the latter over time, with C„ kept as 
quasi-integrals, we get 

E \Cnit)\' = 1 . (21) 
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The quantity 

p„(t) = |C„(t)|2 

plays the role of the mode probability, or fractional mode population, which is normalized 
to one, according to Eq. (21). 

Substituting expansion (19) into Eq. (4), we use the averaging techniques [52-54]. The 
resonant field (13) can be written as 



V{r,t) = - B{r)e"^' + B*{r)e 



'iujt 



where 



We need the notation for the interaction amplitude 

otmn = ^^o'^o / lv^m(r)|^ 2|v9„(r)|^ - |(/9m(r 



dr 



(23) 
(24) 

(25) 



pumping-field amplitude 

I3mn = Jv*mir)Bir)ipr,ir)dr, (26) 

where B{r) is given by Eq. (24), and for the quantity 

7,„ = a,„ - <l>oy'<(r) |2 [pi(r) - pi(r,t) ] ^„(r) + (7i(rX(r) + i^ Ur . (27) 

Also, we introduce the effective detuning 

^mn ^ l^ CO ~r Oimm Oinn ; \^^ ) 

in which 



ann = No^o / lv^n(r)|'^ dr . 



From expression (27) it follows that there exists an effective time te//, during which 7„„ 
can be treated as a real quantity, such that 

|Im7„„| te//<l, (29) 

being weakly dependent on time, in the sense that 



inn (^'^ 



This effective time is of the order 



^e// 



< 1 . (30) 



(31) 



where pi is the density of uncondensed atoms. When practically all atoms are condensed, 
so that pi ^ 0, then tg// -^ oo. Defining 

Cn(t) = Cn(t) exp(i7„„t) , (32) 

we see that, for the times shorter than te//, the fractional mode population (22) can be 
written as 

Pn(t) = |c„(t)|2 (0<t<te//). (33) 

The initial condition (11) for the mode amplitude (32) takes the form 

c„(0) = 5„o . (34) 

With this initial condition, we obtain the equations 

dt 2 

t ^ = a,o\co\'c, + 1 /3,Voe-^^-* , (35) 

where Aqi = A c<j + aoo ~ ciii- 



Equations (35) can be simplified by separating the absolute values |c„| and the phases 
7r„ of the complex quantities 

c„ = \cn\e'^-' . (36) 

Let us also introduce the following parameters 

1 _• 1 

a = 2 ("01 + ttio) , /5 = l/^oil = /^oie '^ , 6 = Aqi + -(aoi - "lo) • (37) 

By defining the population imbalance 

I |2 I |2 /oo\ 

s = \ci\ — \co\ (38) 

and the phase difference 

X = TTi - TTo + 7 + Aoi , (39) 

we can transform Eqs. (35) to the two-dimensional dynamical system 

ds n n . dx I3s ^ , ^^. 

— p V 1 — s^ smx , -— = as A — . cos a; + o . (40) 



dt dt y/T 



s 
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Solving these equations defines the fractional mode populations (33) as 

Po{t) = -^ , Plit) = -^ . 

Equations (35) and (40) were derived earlier [15-17] for a purely coherent system at zero 
temperature, when all atoms were in BEC, so that A^o = N. Here we have showed that the 
same equations can be obtained for a system at finite temperature, when Nq < N. The main 
difference is that, when the density of uncondensed atoms is not zero, then Eqs. (35) and 
(40) are valid not for all times, but in the time interval < t < te//, limited by the effective 
critical time (31). 

In the same way, we could derive the equations for the dynamics of several topological 
modes, generated by the quasiperiodic modulating field 



V{r,t) = W[Bnir)e''"' + B: 



'r)e-ie^t 



with several frequencies e„ tuned to the resonance with different transition frequencies 
Umn = Em — En- FoT example, the equations for three topological modes would have the 
form as in Refs. [31,32], or similar to the equations for three coupled BEC [55]. The systems 
with multiple generated topological modes display a variety of interesting effects, such as 
interference patterns and interference currents [20,24], mode locking [15,24,26], dynamical 
transitions and critical phenomena [17,20,21,24], chaotic motion [31,32], harmonic gener- 
ation and parametric conversion [31,32] that are analogous to these effects in optics and 
for elementary excitations in Bose-condensed systems [56-58], atomic squeezing [24,27,28], 
which can also be called spin squeezing, Ramsey fringes [59], and massive entanglement 
production, which, being similar to the entanglement of two atoms, differs from the latter 
by occurring for multiatomic condensates. In the following sections, we explain how it is 
possible to create and regulate entanglement in a Bose-condensed system with topological 
modes in optical lattices. 



4 Coherent states 

First of all, we need to define a basis of states that we shall consider in what follows. It is 
convenient to use the basis of coherent states. 

Let us consider an optical lattice with Nl sites, each site representing a deep well with 
a large filling factor Vj ^ 1 and with the number of condensed atoms in a well AT,-, so that 

N, = Y.N,, N = J2u,, (41) 

j j 

where j = 1, 2, . . . , N^. Suppose that each well is subject to the action of a modulating field 
Vj{r,t) generating topological modes inside that well. Let a multi-index Uj label the topo- 
logical modes in the j-th site well, and let ?7n(r) be the related coherent modes, normalized 
to the number of condensed atoms in the well, 

N,=J\vnM'dr- (42) 

For generality, we consider the case when Nj and Uj can be different for differing lattice sites. 
This can happen, e.g., if the lattice is perturbed by a disordering potential. 

Topological modes are the solutions to nonlinear equations of type (9), because of which 
they are not necessarily orthogonal to each other, so that the scalar product 

N.J ^ I VU^) Vn.ir) dr (43) 

is, generally, not zero for i ^ j. The diagonal elements of Njj = Nj are the numbers of 
condensed atoms (42). 

For an Uj-mode, we may construct the coherent states in the Fock space as 



\nj > 



exp{-Nj/2) 



k 



UVnM) 



1=1 



(44) 



which is a column with respect to k = 0,1,2, .. .. The coherent states (44) are not orthogonal 
to each other, yielding the scalar product 

< ni\nj > = exp '-- — ^ + A^'^^ , (45) 



but each of them is normalized to one, so that < nj\nj >= 1. 
Let us define the correlation factor 

A.. = ^^ , (46) 

for which Xjj = 1. From the Cauchy-Schwartz inequality for Eq. (43), we have 

\N,,\'<N,Nj (t^j). (47) 

Hence, for factor (46), we get 

8 



The expression in the exponential of the right-hand side of Eq. (45) can be rewritten by 
using the equahty 

liN, + N,)^N,^\(N.-N,f^(l-X,)N.N,. 

The latter diverges if either Ni or A'^-, or both, tend to infinity and i ^ j. Thus, we come to 
the conclusion that the coherent states (44) are asymptotically orthogonal, 

< mlrij > ~ 6ij {Ni + Nj > 1) . (48) 

They also are asymptotically complete (or overcomplete) in the weak sense, 

Y, \nj >< nj\-l {Ni + Nj > 1) . (49) 

rij 

Therefore, the set {\nj >} forms a basis, which is asymptotically orthogonal and complete. 
The closed linear envelope of this basis forms a Hilbert space Tij. The tensor products 

|n > = ^jlrij > (n = {rij}) (50) 

compose an asymptotically orthogonal and complete basis {|n >}, whose closed linear en- 
velope is the Hilbert space Ti = ^j'Hj. The states of BEC, which is a coherent subsystem 
of the physical system, can be interpreted as the vectors of the space Ti. 

5 Lattice register 

By varying the resonant modulating field acting on the lattice, it is feasible to govern the 
creation and behavior of the topological coherent modes and to regulate entanglement pro- 
duced in the system. To quantify the level of the produced entanglement, we shall use the 
measure of entanglement production introduced in Ref. [60]. 

The density operator, characterizing the coherent modes, can be represented as an ex- 
pansion over the basis {|n >}, 

P = Y^ Pn|n >< n| , (51) 

n 

with the normalization 

n 

Let us define a single-partite operator 

Pj = TrH\w, P ■ (52) 

From this definition and Eq. (51), we have 

Pj = Y Pn l^j ><nj\ . (53) 

n 

Then we define the factor operator 

P® = ®fp, , (54) 



for which 

Tr^/)^ = nTrH,P, = 1. 
j 

The measure of entanglement, produced by the density operator (51) is defined [60] as 

e{p) ^ log ^^ , (55) 

\\P \\d 

where the logarithm is to the base 2 and || ■ H^i implies the norm over the disentangled set 

D = {f = ®,v^, |v9, G n,] . (56) 

Here the norms are defined as follows. The set D is assumed to be unitary, that is, for 
any two vectors f & D and /' G -D, one can introduce the scalar product (/,/'), which 
is the standard requirement for any physical system. Having the scalar product makes it 
straightforward to define the vector norm 



generated by this scalar product. Then, for any linear operator A on Z^, the operator norm 
is given as 

\\a\\d = sup , III': "1/21 (/^O, /VO). 

/,/' WjWdWj \\d 

This can also be represented as 

\\A\\n ^sup|(/,i/0| (||/lb=||/'||D = l). 

Vectors f (z D and f E D have the product form similar to Eq. (50). 
For the norms of operators (51), (53), and (54), we have 

\\p\\D = snppn, WpjWn, = snp ^ Pn , ||p®lb = H HPillw. • 

"^ n(^n,.) j 

So that for measure (55), we obtain 

^(P) = log T^ ^^ • (57) 

Entanglement is generated in the system if and only if 

suppnT^n^^^P Yl Pn • (58) 

j "^ n(^n,) 

For example, if all lattice sites would be completely independent, such that p^ would be a 
product of some Pn^ , Pn — ^ Yij Pn, ■ Then, since 

supHPni -^ Ylsnppn^ , 
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measure (57) would be zero, e{f)) -^ 0, that is, no entanglement would be produced. 
The opposite case would be if all lattice sites were correlated, so that 

Pn = Pn n ^""j • (^9) 

3 

This can be realized if all the lattice is shaken synchronically, with the same topological 
mode being generated in all lattice sites. In that case, 

SUpPn = p„ , XI Pn = Pn^nnj , 



because of which 



I ;^ I I — Ti I I ^® I I — r,^L 

\P\\d — Pn , \\P \\d — Pn 



As a result, the measure of entanglement production (57) becomes 

e{p) = (1 - Nl) log sup Pn ■ (60) 

n 

If M topological modes are simultaneously generated in each lattice site, then the maximal 
entanglement production is achieved for p„ = 1/M. In such a case, keeping in mind that 
Nl ^ 1, one gets 

€{p) = Nl log M . 

For the two-mode case, this reduces to e{p) = Nl log 2. By varying the resonant modulating 
fields, it is possible to regulate entanglement in a wide diapason between zero and Nl log M. 
The fractional mode populations are given by p„ = Pn(t), which is defined in Eq. (33). 

The resonant process of mode generation is a fast process, occurring on the time scale 1/a. 
The latter is much shorter than the thermal effective time (31), provided that the number of 
condensed atoms Nq ^ A^i is essentially larger than the number of uncondensed atoms A^^i. 
Another temporal restriction is imposed by the power broadening, defining the resonance 
time tres, after which nonresonant levels become excited, even though the modulating field 
is resonant. The resonance time can be estimated [24] as 






a'^uj 



/32(a2+^2) • 



For P < a and a <C u;, the resonance time tres ~ ^ I P"^^ is much longer that 1/a. It looks, 
therefore, feasible to achieve sufficiently long decoherence times allowing for the functioning 
of the lattice register that can be used for quantum information processing and the creation 
of a boson lattice quantum computer. 

In conclusion, we have shown that the generation of topological coherent modes is feasible 
in Bose systems not only at zero temperature and under asymptotically weak interactions, 
when the whole system would be almost completely condensed, but also at finite temper- 
atures and interactions. This becomes possible because of the resonant character of the 
suggested mode generation. An important feature of the resonant mode generation is the 
feasibility of controlling the process, thus, allowing one to govern the level of entanglement 
production realized in an optical lattice. Such a possibility of regulating entanglement pro- 
duction in a lattice could be employed for creating boson lattice registeres for quantum 
information processing. 
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